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Using the relativistic version of Landau theory of Fermi liquid with a — uj and p mesons exchange, 
we have obtained an equation state for dense neutron star matter in presence of strong quantizing 
magnetic field. It is found that in this scenario the self-energies of both neutron and proton compo- 
nents of dense neutron star matter become complex under certain physical conditions. To be more 
specific, it is observed that in the exchange diagrams of a, uj and neutral p transfer processes and in 
• the direct interaction diagram with p± transfer reactions, the nucleon self-energies become complex 

f"^ ' in nature. 
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1. INTRODUCTION 

^: 

■ With the recent discoveries of a few magnetars 0, E El H- the study of dense neutron star matter in presence 
t— I | of strong quantizing magnetic field has gained a new life. These strange stellar objects are supposed to be strongly 
magnetized, relatively young (age ~ 10 4 yrs.) neutron stars and are also believed to be the possible sources of 
anomalous X-ray pulses and soft gamma ray emissions (AXP and SGR). As surface magnetic fields of these high 
magnetically-powered neutron stars are observed to be > 10 14 — 10 15 G, then a theoretical estimation from virial 
' theorem predicts that the field at the core region may be ~ 10 18 G. Now, if the magnetic fields are really so strong, 
ly-j , particularly at the central region, they must affect most of the physical properties of those strange stellar objects. 
Q\ ' Also, different physical processes taking place at the core region should get modified significantly. 

The strong magnetic field at the core region affects the equation of state of the neutron star matter and as a 
in . consequence, some of the gross properties, e.g. mass-radius relation, moment of inertia, rotational frequency etc. of 
neutron stars should change to a great extent HIS 0> HI • 

In the case of compact neutron stars, the possibility of first order phase transition to quark matter at the core 
Qh region is found to be completely forbidden if the magnetic field strength exceeds 10 15 G. However, a second order 
q , phase transition is allowed, provided the field strength is < 10 20 G p^Hof. 

!_i ■ We have noticed that the elementary processes, e.g. weak and electromagnetic reactions and decays are also 
"o^ \ strongly influenced by such intense magnetic field. As the cooling of neutron stars are mainly controlled by neutrino 
emissions through weak processes (URCA or modified URCA), the presence of strong quantizing magnetic field 
should significantly affect the thermal evolution of neutron stars. Further, the rates of the electromagnetic processes 
• ""H , in presence of strong quantizing magnetic field should modify significantly the transport pro perties of dense stellar 
matter, in particular for electron gas present at the core region of the neutron stars [TTL Il2l ITU . IbH Il5j . We have 
noticed that the electron gas behaves like a neutral super fluid if the magnetic field strength exceeds 10 16 G [l6|. 

It has also been observed that chiral symmetry gets violated in presence of strong magnetic field. The strong 
quantizing magnetic field acts like a cata lyst to gen erate the mass dynamically. A lot of investigations have already 
been done in this field HI El ES M, mil Hill |||. 

In the general relativistic model some investigations have also been done to study the stability of neutron stars in 
presence of strong magnetic field. It is found that there are possibilities of geometrical deformations from the usual 
spherical shapes and in the extreme case of ultra strong magnetic field, they become either black disks or black strings 

In one of our recent studies, we have developed a relativistic formalism of Landau theory of Fermi liquid for dense 
neutron star matter in presence of strong quantizing magnetic field with a — uj meson exchange |29j . In the present 
article we shall re-formulate our previous work as cited above, incorporating p-meson (both charged and neutral 
components) exchange interactions. To obtain the sclf-encrgies for both neutron and proton components from the 
Landau quasi-particle interaction functions considering nuclconic interactions via the exchange of p± mesons, we have 
also included the effect of strong quantizing magnetic field on charged meson propagators. 

To evaluate the self-energy of Landau quasi-particles, we have considered a typical many body fermionic system (in 
this case it is neutron star matter) in presence of strong quantizing magnetic field at zero temperature, interacting 
via scalar, vector and iso- vector meson fields (as mentioned above). This formalism is also useful to obtain the quasi- 
particle energy of other fermionic many body system, e.g., electron gas or even for exotic quark matter in presence 
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of strong magnetic field. The formalism is also applicable to study the magnetic properties and the non-equilibrium 
properties (e.g., the transport properties) of dense stellar plasma of fcrmions. Now similar to the conventional non- 
relativistic model of Landau theory of Fermi liquid, this modified version also deals with normal Fermi liquid and is 
applicable for the low- lying excited states of the system, which consists of a superposition of quasi-particlc excitation 
close to the Fermi surface. The formalism we have developed in this article is basically a re-formulation of our previous 
work 29] and an extension of a very old paper by Baym and Chin 30] . Of course, in the present investigation, we 
have found a new and interesting result, the complex nature of nucleonic self energies in the mean field approximation. 

The paper is arranged in the following manner. In the next section we are giving a general overview of the behavior 
of a fermionic system in a strong quantizing magnetic field. In Section 2, we shall consider a typical dense neutron 
star matter in presence of strong quantizing magnetic field interacting via a — lo — p mesons exchange and develop 
a formalism to obtain the equation of state for dense neutron star matter in presence of strong quantizing magnetic 
fields. In the case of p± and p2 (or p±, the positive and negative components respectively) exchange interactions, we 
shall incorporate the effect of strong quantizing magnetic field on their propagators. In Section 3, we have concluded 
the results and discussed the future perspective of this work. 



2. BASIC FORMALISM 



1. Charged Fermion in Strong Magnetic Field 

In the case of a typical many body fermionic system, let us consider the ith type fermion of charge qi and mass 
(of course, in the case of neutron component of neutron star matter, since qi = 0, then one has to use the original 
work by Baym and Chin and the same formalism is also applicable if the magnetic field strength is less than the 
quantum critical value). We assume the gauge = (0, 0, xB m , 0), so that the constant magnetic field B TO is along 
z-axis. If the magnetic field strength B m exceeds the quantum critical value i?^ = mf/qi, the Landau levels for the 
ith charged component will be populated and the matter is called the Landau diamagnetic system. In the relativistic 
region the charged matter shows Landau diamagnetic behavior, if the cyclotron quantum of the charged components 
exceeds the corresponding rest mass energy or equivalently the de Broglie wave length exceeds the Larmor radius of 
the charged particle. Under such circumstances the phase space integral will get modified and is given by 



fj, 



(2?r)3 j d* P F( P ) = ^ J d Pz d^F{ P ) 



*^£(2-f,J / d Pz F{ Pz ,v) (1) 



where gi is the degeneracy, v = 0, 1, 2, ... is the Landau quantum number and F is a function of particle momentum. 
At zero temperature, the above modified form reduces to 
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i^ l g(2-^ ) f +PFi d Pz F(p z ,v) (2) 



where 



WSL] - 0) 

is an integer but less than the actual value of right hand side, p,i is the chemical potential of the ith charged component 
and P Fi is the corresponding Fermi momentum. Now it can very easily be shown from the solution of Dirac equation 
in presence of such strong magnetic field that the zeroth Landau level is singly degenerate, whereas all other levels 
are doubly degenerate (which is incorporated in eqn.(l) by introducing the factor 2 — S^o)- Further, the energy eigen 
value is given by 

E v = (p 2 z + m 2 + 2vq l B m ) 1 ' 2 = (p 2 z + m 2 + P \^' 2 (4) 
and therefore the chemical potential is simply given by 

Mi = (Pf-i +m 2 + 2vq i B m ) 1 l 2 (5) 



3 



with pj__i = {2vqiB m ) 1 / 2 the orthogonal component of charged particle momentum (in the x — y plane), which varies 
in a discrete manner, hence the name, Landau quantization. 

In presence of strong quantizing magnetic field as discussed above, it can very easily be shown that in standard 
notation, the modified form of Dirac spinors for the charged components are given by 



1 



where 



exp(-iEJ + ip y y + ip z z)u 



U 



(6) 



+m)]V2 



and 



+ m)]V2 



(E v + m)I u:Py (x) 


PzL;p y (x) 

\ -i{2vq i B m ) x l' 1 I v _ Xwv {x) ) 



( \ 

[E v + m)I„_ 1 . Py (x) 
ii2 vqi B m fl 2 I v , Py {x) 
-p z I v -i- p Ax) 



(7) 



(8) 



where the symbols | and [ indicate the up and the down spin states respectively and 



qB n 



1/4 



1 



2 2^/ 2 exp 



\ qB " 



(qB m ) 



1/2 



Py 
_ _Py_ 

qB„. 



(9) 



with H u the well-known Hermite polynomial of order v and L y , L z are length scales along y and z directions 
respectively. 

Now for a dense relativistic many body system of fermions interacting via scalar, vector and iso-vector fields, 
represented by </>, V 1 and p 5 " with masses £, r\ and 9 respectively, the basic Lagrangian density is given by 



1 — -* *~ — — — — T 



+C V 



£>S + £l-V + &EM 



(10) 



where = + iqA^, g v , g s and g p are the coupling constants corresponding to vector, scalar and iso-vector 
interactions; Cy, C-s, C-i-v are the free Lagrangian densities for vector, scalar and iso-vector field respectively and 
Hem is the Lagrangian density electromagnetic part corresponding to constant external magnetic field B m along 
z-axis. Here the iso-vector fields behaves like a four vector in configuration space (represented by yU = 0,1,2,3) but 
like a three vector in iso-spin space (represented by the symbol — >). 

Now in the Landau theory of Fermi liquid, to evaluate various physical quantities of the system, it is a common 
practise to switch on the interactions adiabatically. The normal Fermi liquid system will then evolve from non- 
interacting states to interacting states and there will be a one-to-one correspondence between the dressed or quasi- 
particle states of the interacting system to the bare non interacting states. To obtain /, the Landau Fermi liquid 
interaction function, also known as the quasi-particle interaction function, we have considered both direct and exchange 
type interactions. The exchanged bosons are either scalar or vector or iso-vector type. The quasi-particle interaction 
functions are related to the two-particle forward scattering via the exchange of any one of these bosons. We further 
assume that the external magnetic field acts as a constant background and is present everywhere. In the next 
subsection we consider the processes with scalar boson exchange. 



2. Scalar Boson Exchange 



To compute the Landau-Fermi liquid interaction function / from two particle forward scattering matrix with scalar 
boson (a) exchange, we start with the definition of transition matrix element, given by 



i J ' p s {x)<j){x)d A 



(11) 
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FIG. 1: Basic Elementary Processes 

where 

p s (x) = g s ^(x)^(x) 

is the scalar density and (f>(x) is the scalar field satisfies the the well known Klcin-Gordan equation 

(d 2 + m 2 )<j)(x) =p s (x) 

The solution of this equation can be obtained by Green's function technique and is given by 

d 4 q exp[— iq.(x — x') 



-p s (x')d 4 x' 



q 2_ ( 2 



(12) 

(13) 
(14) 



where q^ = (q°,q) is the exchanged four momentum. 

Now to evaluate Tfi, we substitute the Green's function solution (cqn.14) into eqn.(12) to obtain p s (x) and finally 
putting it into eqn.(ll). We consider the following strong interaction processes: 



(i) p + p- 

(ii) n + n 
(Hi) p + n - 



p + p, 

» n + n and 
> p + n. 



The basic diagrams are as follows: We evaluate Tfi for both direct and exchange type interactions for all these 
processes. We first consider the scattering process represented by (i). The transition matrix element is given by 



_ f , . d q expf— iq(x — x')] . .. jA jA , 
Tfi = ~i Ps(x) 7KZU Pl J y „ ^p^x'^xd^x' 



e 



(15) 



Now using the up and down spin positive energy spinors for proton, we have (since the temperature of the system -C 
the typical chemical potential for nucleons ~ lGeV, we have not considered the negative energy states). 



Ps{x') = — ex$[-i{{E Vl - E V2 )t' - (jpiy - k ly )y' - (p lz - k lz )z'}][p s (x')} 

Li,, Ll 



(16) 
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Then substituting for (f>(x), p s (x) and p s (x') and integrating over t', y' and z', we have the direct part for p — p 
(represented by («)) scattering matrix element. 



fi,(pP^PP) 



= -i f i^8(q - E Ul + E„ 2 )5(q y - p ly + k ly )5(q z - p lz + k lz ) 



exp{iq.x) 
q 2 -C 2 



[p( x )] T 9>> T exp(iq x x')dx' 9s [p(x)]rf 4 a: 



exp[z{(£^ - E„> 2 )t - (p 2y - k 2y )y - {p 2z - k 2z )z}] 
Next integrating over t, y, z and q a , q y , q z we have 



(17) 



rp(d,s) 

fi,(pp^pp) 



-i(27r) 3 S(E)6(p y )6( Pz 



-dxdx' 



(L y L z r 



^[p(x)][p(x')\ exp[-iq x (x - x')] 



q 2 -C 2 



(18) 



where the S functions represent the abbreviated form of energy and momentum conservations along y and z directions. 
Finally, evaluating the contour integral over q x and using the conditions for forward scattering given by E Ul = E V2 = 
E v (say); E^ = E v > 2 = E v > (say); p ly = k ly = p y (say); p 2y = k 2y = p' y (say); p lz = k lz = p z (say) and 
Viz = k 2z = p' z (say), we have the Landau- Fermi liquid interaction function corresponding to direct interaction 



f(d,s 



^ J cxp[-C\x - x'\}[p(x)}[p(x')}dxdx' 



(19) 



where 

[p s {x)] = u(x, u 2 , ki)u(x, vi,pi) (20) 
Then summing over the final spin states and averaging over the initial spin states we have for direct diagram 



\Ps{x)] = ^Tr[A(o;,i/i,i/2,pi,fci)] 



(21) 



and similarly for x' . Whereas for the exchange diagram 



[p s (x)][p s {x')\ = -Tr[A(x,i> 1 ,V2,Pi,ki)A(x, v[ 1 v 2 i' ,p 2l k 2 ) 



(22) 



where 



A(x,v 1 ,v 2 ,p 1 ,ki) = y~] u(x, vi,pi)u(x, v 2 ,k 2 ) 



(23) 



To obtain A's for the process indicated by (i), we consider the direct diagram for forward scattering. Then 
substituting the up and down spin solutions of Dirac equation for proton (eqns.(7) and (8)) in the expression for 
projection operator (eqn.(23)), we get 



A = —^—[Ak^-f^ip = and z) + mA + Bk^(p = y and p y = p ± )} 
Lb,, 



(24) 



The matrices A and B are given by 



(I V I' V 

I v -il' v _ x 

o o i v r v o 

Vo o I v -xl' v _ x 



(25) 



B 



( i.-iK 











T J' 











h-1 


Vo 














(26) 
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where the primes indicate the functions of x' . To the best of our knowledge, eqn.(24) is an entirely new result, first 
obtained in our previous publication 29] and has not been reported earlier. Further, these results play key roles in all 
kinds of calculation related to the electromagnetic and weak interactions in presence of strong quantizing magnetic 
field. In a recently completed work 25] we have used this results to study chiral symmetry violation in presence of 
strong magnetic field. 

Since 7 matrices are traceless and both A and B matrices are diagonal with identical blocks, it is very easy to 
evaluate the traces of the product of 7-matrices multiplied with any number of A and/or B, from any side with any 
order e.g., 

Ti(yi J -j v A 1 A 2 ..B 1 B 2 ..) = Ti(A 1 A 2 ..B 1 B 2 ..)g>"', (27) 



Tr(yV 7 VAA2..Bi£ 2 ..) - Tr(A 1 A 2 ..B 1 B 2 ..)(g^g- x - g^g va + g^g va ), (28) 

Tr(product of odd 7s with A and/or B) = etc. The other interesting aspects of A and B matrices are: 

i) k lfl k 2 ^Tr(A 1 A 2 ) = (E 1 E 2 - k lz k 2z )Tv{A 1 A 2 ) 

ii) kx^TtiBxBi) = k 1± .k 2± Ti(B 1 B 2 ) 

iii) fc^A^Tr^!^) = kin&TriBtAi) = 

iv) p^k^p^TriA^) ± = {E V1 E V , 2 - pi^p^M^A^) 

These set of relations are also totally new results and have been derived in our previous publication |29| . 
Then we have following Baym and Chin [30| 

f£%p = - 16 f£ Ev , J ex P[~C \ x ~ x'\]TrATvA'dxdx' (29) 

where 

TrA = 2{ll + Il_ x ) and TrA' = 2(# + ( 30 ) 
On the other hand if we consider the direct interaction corresponding to the process (ii), we have 

f{d,s) = 9s m2 /o-|l 
Jnn—>nn C, 2 EE' 

This result is same as that obtained by Baym and Chin long ago. 

Next we consider a bit complicated case; scattering process represented by (iii). In this case only protons are 
affected by strong quantizing magnetic field, their Landau levels are populated. The transition matrix element is then 
given by 

T fi = i J pP(x)^(x)d 4 x (32) 

where pi P \x) is the scalar density corresponding to proton component and is given by eqn.(20). The scalar field 

**) = " / ^ P[ ~^ X ' )] pi n) (^ (33) 
where the scalar density corresponding to neutron component is given by 



where 



pf\x) = g s i> n (x)i> n (x) = ^exp[-z( P2 - k 2 ).x][pf\ (34) 



[pW]=u(k 2 )u(p 2 ), (35) 
here u and u are the usual form of spinor and the adjoint for neutron, given by 



u 



A£j V E+mX. 



(36) 
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where s = 1 and 2 corresponding to the up and down spin cases respectively, and are given by 

X (1) = X T = ( M (37) 



and 



X (2) = X X = ( J ) (38) 
Using all these expressions and integrating over x, we have 

^• ( i_ p+ „) = -*/ rf^^^l^l^exph^^ - E V2 )t' i(p ly k ly )y> i(p lz - k lz )z'}\ 

[p ( i' ) ]5 A (q-p2 + k 2 )d 4 x' (39) 
Further integrating over t' , y' , z' and q x , q y , q z , q finally we have for the forward scattering case 

fp+n^p+n = - ^2^J^ m ( 4 ^) 

where we have used the same expression as before for [p^ (x)] whereas 

lp { s n) ] = 7^(P^ + m) (41) 
Here we have also used some obvious results as given below to obtain the above expression (eqn.(40)). 

/+oo 
TrA'dx' = 4 (42) 
-oo 

Although the structural forms of both eqns.(40) and (31) are identical, the appearance of proton energy in the 
denominator of eqn.(40) makes it qualitatively different from eqn.(31). This is also obvious from the physics point 
of view. The interaction with c-meson exchange can not distinguish neutron and proton. The difference comes from 
the change in proton parameters (spinor wave functions, energy eigen values etc.) in presence of strong quantizing 
magnetic field. As we will see later that in the evaluation of energy densities for proton and neutron separately the 
form of momentum integrals (phase space integrals) will differ significantly (which has also been discussed in the 
previous section). 

We shall next evaluate the Landau interaction function for exchange part with cr-meson exchange. In this case the 
processes represented by (i) and (n) give results in which there are nothing surprising. Whereas, as we shall show that 
the Landau interaction function obtained from the process (Hi) is complex in nature. Which will give complex self 
energy of both proton and neutron. This is also true for exchange diagrams with vector boson and neutral iso-vector 
boson exchange, whereas in the case of charged iso-vector exchange, the direct term gives complex Landau interaction 
function and as a result the complex self-energies. The findings are to the best of our knowledge are new and not 
been reported earlier. 

We start with the process (i). In this case the forward scattering condition gives E Vl = E v ^ = E v (say), E v ^ = 
E V2 = E v , (say), p ly = k 2y = p y (say), p lz = k 2z = p z (say), p 2y = k ly = p' y (say), p 2z = k lz = p' z (say). 
Then we have the Landau interaction function 

= / dxdx' eM - k ^- A] [{E v E v ,-p z p zl+ rn^{AA') 

-pl.p 7 ± Tr(BB')} (43) 

where K = (q y + q 2 z + C 2 ) 1 ^ 2 , q% = {p% — p'j) with i — y and z. The traces are obtained from eqns.(25) and (26). 
Similarly for the process (ii), the Landau quasi-particle interaction function is given by 

Aex,a) g 2 p-p' + m 2 

Jn+n^n+n ^EE' (j)' - p) 2 - C, 2 K ' 



and is identical with zero field result of Baym and Chin. 
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Next we consider the exchange term corresponding to the process (Hi) with scalar boson transfer. In this case the 
transition matrix element is given by 



T/^y p<r>(x)<f>(x)#x 



(45) 



where 



(x) 



3 ip p (x)ip n (x) 



= t~T exp[-i{(-B P - E k )t - (p y - k y )y - (p z - k z )z}} exp(ip x x)[p { f n \x)} 



LyL z 



and 



Hence 



m = - 



d 4 q exp[— iq(x — x')\ 



(2n) 4 



q 2 -C 2 



p^ p \x'd 4 x' 



(46) 



(47) 



T fi = -i 



d q 



(2ir) 4 L y L z 
cxp[-iq(x - x')] g. 



exp[-i{(2^ - E v ,)t (p ly k 2y )y (p lz - k 2z )zWr\x)] 



LyL z 



exp[i{(E^ - E V2 )t' - [p 2y - k ly )y - (p 2z - fcijz}] exp[-i(pi x - fci Jx]d 4 xd 4 x' 



(48) 



Now integrating over t, y, z, t' , y', z', qo, q Xl q v and q Zl we have the Landau quasi-particle interaction function with 
the condition of forward scattering 



(ex.s) 



4/ 



2 f ^ exp[ ^ g --^^ + /-l ( .f- a;/)] TV[A p „(^P)A, P ^^0] 



Ip+n^p+n - 4 J — ( p _p')2_^2 

Now to evaluate the projection operators A np and A p „, we define 

Px + ip y = P± exp(i<f>) and p x - ip y = p± exp(-i(j>) 



(49) 



(50) 



Then from the ranges of p x and p y , both of which vary from — oo to +oo and the range of p±, which is given by 
(2vqB m ) 1 / 2 , with < v < ^ max where f max is a function of density of matter and the strength of magnetic field, we 
have < < ir. With these choices, we have 



K P (x,p) = ulul + u^up 

( I v {E + m) 

1 

2E 





Vzh 



I v -i(E + m) 
p±e- i +I v - 1 



-Pzh 

-ip±L 



-ipi_I v -\ \ 
PzL-i 



(ipje-^-pj) T PzP±(e- i < l >-i) T 

lv c 1 — +V — 1 



E+m 



E+m 



(51) 



\P±d lv Pz^v-l E+m iv E+m J-v-l / 

Which after some simplification reduces to 



Anp " 2E 



(Ej + m-p z ~f z )A- 



- -S z (l -7o) - ^ ; — (l-7o)S„ z ^. 



£ + TO 2 



2(£ + to) 



" -m 2 k 
(1 + 7o) 



} j x A-p±-y x — {l + ~/o)X* 



where 



X 



( e-** \ 



e"^ 

V J 



(52) 



(53) 
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(JiOj U 
<7j<7,- 



(54) 



I v 0\ 

I v -y 

J„ 

I v -xJ 



(55) 



and 

<Ji 
^ 

whereas, 7's are Dirac matrices in standard representation and <Tj are the usual Pauli spin matrices. Similarly 



(56) 



1 

2E 



( I v (E + m) 
h-xiE + m) 



-Pzh 
-p^h-x 

[(E 2 -rn 2 )-pl(te^-l)] 



-P±e-^I V \ 
PzL-i 



E+m 

2 



(57) 



E+m 



j bl(»e-^+l)-(^-m^)] j I 

±v ~ 1 E+m ±v ~ 1 / 



After some simplification as before this can also be written as 



A J_ 

pn ~ 2E 



2 

(£70 + m - Pzlz )A + ./| , (1 - 7o)£ z »4 + , (1 - l0 )V x A 



+ 



2(E + m) 



2(£ + to) 



£ + to 2 



(1 - 7o )X* 



E TO X — (1 - 70) 
£ + to y 2 1 



+ pi_AH x X* - (1 , 7o) S^ 



(58) 



Since diagonal elements for X^, E,j and 7, with i ^ j and i,j = are zero, the traces of all the quantities 

multiplied with these matrices are identically zero. Therefore in both A np and A pn only the first and the fourth terms 
will contribute in trace calculation. As a result the Landau quasi-particle interaction function will become complex 
in nature. The explicit form of the trace term is given by 



Tr [A np (p')A pn (p)} = ——{(E V E V , +p z p' z + m 2 )2{I v I v > + J„-i/„'-i) 

A{I V I V , + J„_iJ„,_i) 

iEvp'l 



2 / 

p±p± 



4(E v + m)(E v , +m ) 



2{E„. + m) 



imp 



2{E V > + to) 

iE v >p\ 
2{E V + to) 
I v -iL'-i exp(z</>))] 



2{I V I V , exp(i^) + 7„_i J„,_i exp(-i^)) + 



2(l v I vl exp(i(j)) + I v _\I v i_\ exp(-i(f>)) 



2{I V I V , exp(-i</>) + /„_!/„,_! exp(z0)) + 



ITTiPj 



2{E V + to) 



2{I V I V , exp(-i<j>) + 



(59) 



Obviously the last four terms in the above expression are complex in nature. 

Now the self-enegies for proton corresponding to the direct and exchange diagrams of process (i) is given by 

[ ^max ] 

4" /ei) fc) - E v {p z ) + £ (2 - <W) / / d Pv d Pv ,dp x ,fW%£) 

{ I u'=0 Jp' y =-oo Jp z ,=- Pf 



EA Pz ) + AE^ Dii/ J Pz ) 



(60) 
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For the process (ii) the neutron self-energy corresponding to direct or exchange digram is given by 

e W«>(p) = E(p) + ^- 3 J f^ +n d 3 p' 

= E(p) + AE^l /cx (p) (61) 
Next for the process (Hi) with direct diagram, the proton self-energy is given by 



e^(p z ) = EAp z ) + ^f:\2-S^) jf£ S n l P+ ,M 

u'=0 



= E v (p z ) + AE% ix/ J Pz ) (62) 

Similarly for neutron, it is given by 

e^(p) = E(p) + j^y 3 j f^ P+n d 3 P ' 

= E(p) + AE^ r (p) (63) 

Finally the exchange part of self-energy corresponding to process (Hi) are same for both proton and neutron and it 
is given by 

["max] ~ n f+PF 

e^\ Pz ) = E v ( Pz )+ ]T (2-J„ ) / ##W±cos(^)cos(^) / dp'J^l p+n 

= E v (p z ) + AE^(p z ) (64) 

This last expression for self energy for both proton and neutron is complex in nature. As we will see in the subsequent 
sections that similar results can also be obtained with vector and iso- vector exchange type interactions. 

3. Vector Boson Exchange 

We shall now consider the vector boson exchange type interaction. The transition matrix element is given by 

T f i = iJj»(x)V li (x)d! l x (65) 

where 

f(x)=g v i ) (x) 1 ^(x) (66) 
is the particle current and the vector V^(x) satisfies the Klein-Gordan equation given by 

(d 2 + r ] 2 )V f "(x) =f(x) (67) 
Hence the Green's function solution for V fl (x) is given by 



r) 4 q 2 - rj 

where r] is the vector boson mass. 

Let us first consider the direct term of the process (i). Using the spinor solutions for proton (eqns.(7) and (8)), we 
have 

= y^- exp[-i{(£^ - E V2 )t - ( Ply - k ly )y - (p lz - k lz )z}][f] (69) 
where in the direct term we define 

spin 

= lTr[A(x',u u iy 2 ,p 1 ,k 1 )^} (70) 
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whereas in the case of exchange term we use 

[f{x)][j^x')] = ^Tr[A(x,u 1 ,u 2 ,Pi,k 1 )^A(x , ,u[,u 2 ,p 2 ,k 2 )^} (71) 

Similar to scalar meson exchange case as discussed before, we first integrate over t', y', z'; then over go, q y , q z and 
the contour integral is evaluated over q x . The final integrations are over t, y and z. Thus the Landau interaction 
function for p — p (process (z)) forward scattering with vector boson exchange direct term is given by 

2 r- 

ti d +;l p+p = 32 J^ J exp[-v\x-x'\}dxdx'[(E„E v ,-p zPz ,)TrATrA' 

+pI. P 7 l TtBTtB'} (72) 

The above expression for process (ii) with vector boson exchange direct interaction is identical with the zero field 
case and is given by 

&l n+n = ^^j(EE' P.p') (73) 
Next, for the process represented by {Hi), the matrix element is given by 

Tf^z J f [n) {x)V,{x)d 4 x (74) 

In this case the vector field is given by 

V»{x) = - / ^ eXPt T^'V wV (75) 



where 



= 9 ^eM-i(P2-k2).xM n) ] (76) 



and 



J#V) = -fir eM~i{(E Vl - E V2 )t' (p ly k ly )y' ( Plz - k lz )z'}]\j^] (77) 
Further for the direct term 

= l -Tr[u{k 2 )^u{p 2 )} (78) 

and 

= ^[0(fci)7 M «(Pi)] (79) 

Substituting all these terms in Tji and subsequently integrating over x', t', y', z', x, t, y, z and d 4 q, we have the 
direct type Landau-Fermi interaction function for w meson exchange as 

„2 r + ao 

f^:i P+n = iE X^ p - p ' J x dx ' Tr{A + B) m 



where we have used 

[i(n)W] = ^ Tr [(i' + m)7' t ] = ^ (81) 
and 

[j#V)] = J^ Tr [(V7A + Am + Bp x 7x )^\ = -^Tr{A + B]p^ (82) 
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It can now be shown easily that 

f+oc c + oo 



dx'TrA(x') = 4, / dx'rmTrB(x') = 

-oo J — oo 



(83) 



Therefore 

/^:l p+ „ = -^^(E'E V - (84) 

The structural form is same as that of non-magnetic case but qualitatively the result is entirely different as we have 
already discussed. 

Now we shall consider exchange interaction terms corresponding to w-meson transfer. For the process (i), i.e., 
p + p — > p + p, we have 

j ex ,v) _ al [ + °° [ + °° ^,exj>[-k\x-x'\] 



/p+p ^ p+p " 1GE V V E V , ]_„ ]_„ dxdT K 

[{ PzPz , - E V E V , + 2m 2 )Tr(yL4') - 2p ± .0 ± Tr{BB')] (85) 

In the case of n + n — > n + n, i.e., process (ii) the result is same as that of non-magnetic case and is given by 

Aex.v) al 2m 2 - p.p' 1 

Jn+n^n+n 2 EE , ^ _ pl y _ ^ \°») 

This is exactly the same as that of zero field result. 

Let us now consider the exchange term corresponding to process (Hi). In this case the calculation is almost identical 
with the scalar exchange model. Here unlike the scalar case we have to evaluate Tr[A np (p)j fJ, A pn (p')^^]. where the 
projection operators are given by eqns.(52) and (58). We have noticed that from the non-zero contributions it is very 
easy to show 

Tr[A pn (p) 7 A np (p')7o] - -Tr[A pn (p) 7 a; A„ p (p')7x] - -Tr[A p „(p) 7 ^ 'A np (p') ly ] 
= -Tr[A p „(p) 7 z A np (p')7,] = 

i -2 '2' 

-\{E V E V , + m 2 +p z p' z ) + - J P± <m - E v ,)Tr[AXA'} + — ^ r(m - E v )Tc[AA'X'] 



±E v E v , vy v v 2{E v +m) K v ' r J 2{E v ,+mY 

- 2(E„ + ^, +m) T ^ < 87 > 

where 

Tr[AA'] = 2(I v I v i - I v _\l v t_\) 

Tr[A4'X] = 2(J„J„/ cx P H0) - I v - X I v >-\ exp(i^)) 

Tr[AA'X] = 2{l v l v , exp(i^) - J„_il„/_i exp(-i^)) (88) 

It is quite obvious that complex contribution is also coming to the Landau interaction function and as a consequence 
to the quasi-particle energy if we consider the exchange diagram with vector boson transfer in the process represented 
by (Hi) 



4. Iso- Vector Boson Exchange 

Now, we consider the last part of our investigation; the iso- vector boson exchange in neutron star matter. Since this 
part was not included in the original work of Baym and Chin on conventional rclativistic version of Landau theory of 
Fermi liquid, we shall first develop a formalism for the non-magnetic case. Here, as we have mentioned, the iso- vector 
field behaves like a four vector (indicated by symbol /x) in configuration space and a three vector (indicated by — >) 
in iso-spin space. 
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We start with the simplest case, the exchange of neutral iso- vector boson pg. The corresponding transition matrix 
element is given by 

T fi = if j%(x)p^(x)d 4 x (89) 
where 

*.« - - / ^ M 7 [ S x ' )] ^ l)d,J (90 » 



where 9 is the mass of iso- vector boson and 

= S*)/|W (91) 
is the third component of iso- vector fermion current. In the case of neutral p-meson exchange 

f^) = \\3i (p) ^)-f Un) ^)] (92) 

where symbols p and n represent proton and neutron components respectively. Using the conventional spinor solutions 
for both proton and neutron we have, 

is ( x ) = \ [up(k)Yu p {p) ~ u n (k)^u n (p)} exp[-i(p - k).x] (93) 

On substituting j£ and pg in the transition matrix element and following the procedure by Baym and Chin for 
the case of vector meson exchange, we have the direct type Landau quasi-particle interaction function with neutral 
iso- vector boson exchange 

Ad,iv) = 9 P P-p' , . 

3 46» 2 AEE' [ ' 

Similarly, for exchange type interaction we have 

f (ex,iv) = gp I 2 ™ 2 - P-P'} 

J 8EE> [(p - p' f - 2 ] { ' 

Now in the non-magnetic case, it is very easy to show that the same expressions as above (eqns.(94) and (95)) will 
also be obtained for the charged iso- vector mesons exchange for the valid process represented by (Hi). 

Let us now consider the presence of strong quantizing magnetic field, the Landau interaction functions for neutral 
p-meson exchange are almost identical with the results of vector boson exchange (except a factor of 1/2). For the 
process (z), with direct type interaction, we have 



2 

V 

M6E~E V 

+ P 1. P 7 l TtBTtB'] (96) 



f^i p+p = jg- J cM-e\x-x'\]dxdx'[{E u E vl -p z p zl )^ATrA' 



For the process (ii), the direct interaction is again 1/2 times the non-magnetic result which we have already discussed 
and in the case of (Hi) the Landau quasi-particle interaction function will again be 1/2 of that of vector boson case. 
In the last case, we again got complex contribution in Landau interaction function. 

Next we consider the exchange of charged p-meson. For the sake of convenience we define 

rt{x) = \\fr + {x)+(f-{x)\ ^d p%(x) = ±[ff + (x)-pi*_(x)] (97) 

where p^ (x) are related to the well known Ladder operators in iso-spin space (causes p n) . Then in the lowest order 
diagram, the only valid scattering process is as represented by (Hi). The transition matrix element for non-magnetic 
case is given by 



Tf . - L / ^IL u<Jri( x \ + 4np) (x)) e M-iQ-(x-x')} 

(j? pn) (*)+j? np) ( x » d4xd4x ' ( 98 ) 
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where 



and 



(99) 



(100) 



After some simple algebraic steps we have 



Tfi = -- 



d 4 q 



16 J (2ttY 



[J^ix) 



exp[— iq.(x — x ) 



'j(np)( X ) 



expHg-Qr -x)} 



Using the conventional form of spinor solutions for both proton and neutron, we have 



j(pn)( X ) = y ex P[-*(Pp - Pn)-x]u n (p n )^U p (p p ) 



and 



j^ p \x') = ^ exp[-i(p' n - p' p ).x']u p {p' p )^u n (p' n ) 



(101) 



(102) 



(103) 



Then integrating over d 4 x, d 4 x ; and d 4 g, the Landau-Fermi interaction function for direct type forward scattering 



9pP-P' 
%EE'6 2 



Similarly for the exchange case 



Mv) = 9 2 P ^m 2 +p.p') 
J 16EE'[(p - p'f - 8 2 ] 



(104) 



(105) 



Let us now consider the effect of strong quantizing magnetic field on this process with charged p-meson exchange. 
In this case not only the Dirac spinor solutions of proton, but the charged meson propagators will also get modified. 
The meson propagator is given by 



G p± (x) = 



eB m 
4tt 2 



]T f dE <f>*(x')$(x)exp[-i{E (t-t')-q v (y-y')-q z (z-z')}} 



^- A [jpn{x') - jnp{x')]d A x' 



E - (e q - iS) 4 



(106) 



where 



®(x) = I„{x,p) = 



qB n 



1/4 



exp 



--qB m [ x - 



(qB m ) 



1/2 



Py 

qB m 

Py 

qB„ 



(107) 



Then it is obvious that $*(x) = $(x). The quantity e q = (q 2 z + 6 2 + 2veB m ) 1 / 2 is the modified form of charged 
p-meson energy eigen value, obtained from the solution of Klein-Gordon equation in presence of strong quantizing 
magnetic field, with 2v = 2n + 1 + 2a, a = — 1, or +1 and n = 0, 1, 2, .., the Landau principal quantum number. 

Now in the direct case, unlike the scalar and vector mesons exchange, we have to evaluate the product 
Tr[A p „7 A ']Tr[A„p7 M ]. As before, we calculate this quantity term by term (i.e., for p, = 0,x,y and z). Using the 
expressions for the projection operators, we have 



Tr[A pn 7°]Tr[A„ p7o ] - 



(108) 
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Tr[A p „7lTr[A„ p7 z] = Tr[A p „ 7 ^Tr[A„ p7y ] = 



(109) 



and 



Tr[A pn7 *]Tr[A„ p7z ] 



-^PzTtA 



--p'TrA' 
2 



(110) 



Obviously, the term with /x = is complex in nature and it will make the Landau quasi-particle interaction function 
and as a consequence the quasi-particle energies for both neutron and proton complex in nature. 

Fortunately, in the case of exchange type interaction with charged iso-vector meson exchange, we have to evaluate 
Tr[App 7 M A nn7At ] which gives real contribution to quasi-particle energy. 

To obtain the variation of self energies (both real and imaginary parts) we have evaluated numerically the total 
self energy density for nucleons (both protons and neutrons separately) then plotted the average proton self energy 
as a function of proton density and are shown in the following two figures. To draw these figures, we have solved 
a few non-linear equations sclf-consistcntly. We assume that before the interaction is switched on, the matter is in 
/3-equilibrium, given by 



(Which is obtained from the weak processes n — > p 



and p 



known as the URCA and the 



modified URCA processes). We further assume the charge neutrality condition, given by 



and assuming that the baryon number density tlb 
Hi, where i = n,p and e, we use 



n„ = 



Tip — Tl e 

n p + n n remains constant. To solve for the chemical potentials 

(/4-Q 3/2 

6tt 2 



n p = E^ 2 - & >*M - m l- 2 "^) 1/2 

i/=0 

and 

Hence, knowing the chemical potentials, we obtained numerically the total free energies (both real and imaginary 
parts) per particle for both protons and neutrons. 
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FIG. 2: Variation of the real part of average proton self energy (in MeV.) with proton density (expressed in terms of normal 
nuclear density no = 0.17fm~ 3 ). The upper curve is for B m = 10 14 G, middle one is for B m = 10 16 G and the lower one is for 
B m = 10 18 G. ' 



10 5 = 



10 
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FIG. 3: Variation of the imaginary part of average proton and neutron self energies (in MeV.) with proton density (expressed 
in terms of normal nuclear density no = 0.17fm -3 ). The upper curve is for B m = 10 18 G, middle one is for B m = 10 16 G and 
the lower one is for B m = 10 14 G. 



As one can see from these two figures that the real part always increases with the increase of matter density but 
the matter becomes more stable if the magnetic field is stronger. Which is also true in all other cases, e.g., atoms, 
quark matter etc. |lfl Isil l8^ | in strong quantizing magnetic fields. Whereas the imaginary part decreases with the 
increase of matter density and at very high density the imaginary part is almost negligibly small with respect to the 
corresponding real part. The reason is that at very high density, the effect of Landau levels become insignificant, in 
other ward, at high enough density the magnetic field no longer has quantum mechanical effect. On the other hand 
at low density, it is significant and the effect is more as the magnetic field becomes stronger. 

For the sake of completeness we have also plotted the life time of nucleons (both neutron and proton) as a function 
of proton density, for three different magnetic field strengths. The life times are in the strong interaction time scale 
within a broad density range. Because of such small time scale we belive that the bulk neutron star matter remains 
stable against strong decay, just like a nucleus is a stable object although strong interaction processes are going on 
inside. The other extreme could be that the magnetic field strength at the core region of magnetars are much smaller 
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FIG. 4: Variation of proton and neutron life times (in sec. in unit of 10 22 ) with proton density (expressed in terms of normal 
nuclear density no = 0.17fm -3 ). The upper curve is for B m = 10 18 G, middle one is for _B„ 
B m = 10 14 G. ' 



than the quantum critical value and the observed strong magnetic fields are concentrated at the crustal region. 



3. CONCLUSIONS 



We therefore conclude that in the relativistic version of Landau theory of Fermi liquid in presence of strong quantiz- 
ing magnetic field, although the self-energies for both neutron and proton become complex in nature, the bulk neutron 
star matter still remain stable against strong decay. We have noticed from the complex nature of self-energies for 
both the dressed protons and neutrons, that the life time is within the strong interaction scale and we belive that just 
like a stable nucleus, the bulk neutron star matter will also remain stable in presence of strong quantizing magnetic 
field. 

The other explanation, which is of course not so convincing is that, we can have a stable dense neutron matter at 
the central region of a magnetar with a very low magnetic field strength. The observed strong magnetic field in such 
objects are concentrated at the crustal region. However, in our opinion, the best explanation to have a magnetic field 
of strength ~ 10 14 G at the surface and several orders of magnitude large field at the core is that the conducting (and 
possibly not superconducting because the field strength is > the critical value) neutron star matter does not allow 
the magnetic flux lines to move from the core region to the crust. Since the density of neutron star matter at the core 
region is high enough, a large number of Landau levels will be populated, even if the field strength is extremely high. 
Which further means, that the quantum mechanical effect of strong magnetic field will be extremely weak. 

Finally, we would like to add a few words to compare the results, in particular the complex nature of nuclcon self 
energies we obtained in this work with some of the well known phenomena already discussed in the literature. 

About sixty years ago L.D. Landau had shown the possibility of some kind of damping or dissipation in collisionless 
plasmas and this phenomenon is known as Landau Damping 33, 34]. Being independent of collisions, it is fundamen- 
tally different from dissipation in ordinary absorbing media; the collisionless dissipation does not involve an increase of 
entropy and is therefore a thermodynamically reversible process. In this phenomenon, there is a continuous exchange 
of energy between plasma waves (produced by some external perturbation) and plasma particles. As a consequence 
the self energy of plasma particles become complex and they loose energy to plasma wave reversibly and in the very 
next step again gain energy from the plasma waves and the process of loss and gain continues, but does not decay to 
other particles. This phenomenon makes the system mechanically stable. 

The evaluation of nucleon self energy £ in relativistic approaches has been restricted mainly to the real part of £ in 
either the mean Geld theory or the relativistic Hartree approximation. In these approximations the self-energy is real, 
local and energy independent. It is completely determined once the constant effective mass and the mean constant 
value of vector potential is know j3f| . 

However, the imaginary part of the self-energy, which in reality the imaginary part of optical potential has been 
studied to lowest order by Horowitz |36| . Here, only the on-shell properties of Im(E(e,p)) is important and is 
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restricted to a limited energy range. He obtained a quantity similar to the imaginary part of optical potential by the 
cancellations of large contributions coming from the exchange of scalar and vector mesons; between scalar, time-like 
and three-vector contributions which are also independently very large. 

In ref. real and imaginary parts of nucleon self-energy was obtained with relativistic Dirac-Bruckner approach 
assuming a — u> meson exchange. In this paper fourth order meson-nucleon coupling was considered and studied the 
off-shell behaviour of imaginary part of the nucleon self energy. 

In the reference (3^1 Dyson equation is used to obtain the nucleon self energy £*. In this paper the irreducible 
self energy is calculated by summing over the Dyson expansion to all order. It is shown that in the lowest order 
diagrams the irreducible self-energy contains mean field part and the Hartree-Fock contribution. Both of them are 
real in nature. The dynamic part of the self energy is expressed in terms of the four point vertex function. 

However, the complex nature of nucleon self-energy as we have obtained in this article, arises because of some kind 
of mixing in the iso-spin space (the modified form of Projection Operators e.g., A pn and A np are actually responsible 
for such strange phenomena. These modified form of projection operators are again coming because in presence strong 
quantizing magnetic field, the Landau levels of protons are only populated, whereas neutrons behave classically) and 
it is at present very difficult to predict what will happen to these nucleons with complex energy. 
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